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Abstract. Dark matter particles may be captured by a star and then thermalized in the 
star's core. At the end of its life a massive star collapses suddenly and a compact object 
is formed. The dark matter particles redistribute accordingly. In the inelastic dark matter 
model, an extended dense dark matter mini-halo surrounding the neutron star may be formed. 
Such mini-halos may be common in the Galaxy. The electron/positron flux resulting in the 
annihilation of dark matter particles, however, is unable to give rise to observable signal 
unless a nascent mini-halo is within a distance ~ a few 0.1 pc from the Earth. 
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1 Introduction 

Dark matter (DM) is a form of matter necessary to account for gravitational effects observed 
in very large scale structures such as anomalies in the rotation of galaxies and the gravi- 
tational lensing of light by galaxy clusters that cannot be accounted for by the quantity of 
observed matter [e.g., 1, 2]. The existence of dark matter is a piece of solid evidence for 
new physics that is beyond the standard model and the relevant researches are the focus 
of modern physics and astrophysics. Among the hypothetical particles proposed so far, the 
weakly interacting massive particles (WIMPs) are the leading candidate [1, 2]. The direct- 
detection searches look for the signal of WIMP-nuclei scattering in underground detectors. 
The indirect-detection experiments, instead, aim to catch the annihilation/decay products 
of DM particles. Currently there are some anomalous signals in favor of WIMPs but the 
evidence is not conclusive [3-6]. 

The DM particles may be captured by the stars. The capture of WIMPs by the sun, 
earth and other astrophysical objects has been extensively investigated [7-9]. The captured 
particles will eventually be thermalized and centered in the core of the star, namely, a DM 
core forms in the center of a star [10]. After the death of a (massive) star, a compact object 
is formed. The collapse of the star lasts for a very short time and the density distribution 
of the DM core will be modified. If the compact object is small enough, a significant part 
of the DM core may be left outside the compact object, i.e., a DM mini-halo is formed. 
As the outer shell of the stellar remnant eventually disappears, the annihilation products of 
the DM matter, such as electrons and photons, are visible for the observers. In this work 
we investigate the formation and evolution of such a mini-halo and calculate the flux of the 
annihilated electrons/positrons. This work is structured as follows. In section 2 we describe 
the distribution of WIMPs inside the star. In section 3 we investigate the redistribution of 
WIMPs in the collapse process. In section 4 we calculate the evolution of the DM mini- 
halo and in section 5 we estimate the electron/positron flux detectable on the earth. We 
summarize our results in section 6. 



- 1 - 



2 Dark matter capture by the star 



The WIMPs scatter with nuclei and lose their kinetic energy when pass through a star. These 
particles will get captured if the kinetic energy is smaller than the stellar gravity potential. 
After being captured, the WIMPs will be eventually thermalized in the star [10]. The density 
distribution of the WIMPs is described by [11, 13] 

P{r) = Pcentere-'"-'^^^)/^ = Pc.niore-'" ''"^ (2.1) 



where T is the temperature of the center of the star, (j){r) is the gravitational potential, 
is the rest mass of WIMPs, Pcontcr is the dark matter number density in the center of the star 
and 

27rm^G/)ccntcr l.ikeV Pcentcr 

where Rq is the solar radius. As a result, a DM core is formed and the annihilation efficiency 
is given by [14] 

{Jd^pirW (2vr)3/2.3^' l^-^J 

where < a^v >~ 3 x 10~^^ cm^ s~^ is the averaged annihilation cross section. The number 
of total DM particles (iV) evolves as 

^ = C - CaN\ (2.4) 

where C is the capture rate and the annihilation rate can be expressed as Ta = ^Ca^V^. 
Eq.(2.4) can be solved analytically and we have 

Ta = ^Ctanh2(t/req), (2.5) 

where r,^ = 1/{CCaY/^ = ^^'y^{^?'\^^^S^^f'^{^^^Y'^- ^ ^cq is longer than 
the lifetime of the star, the equilibrium can never be reached. If equilibrium is reached the 
annihilation rate is Ta = i^C. 

In most works, the scatter is assumed to be elastic. However, the recent results from 
DAMA and CDMS may favor the inelastic DM model [15-17]. In such a model, the WIMPs 
will scatter inelastically with the nuclei so the kinetics is modified [18-20]. The cross section 
in the inelastic case can be a few orders of magnitude larger than that in the elastic case. 
Nussinov et al. [11] discussed the inelastic capture of DM particles by the sun and found 
that the annihilation-capture processes can reach equilibrium but the density profile of the 
DM core is very different from that formed in the elastic capture. For example, the radius 
of the DM core is much large. Some DM particles go outside the star [21]. 

3 Dark matter distribution after the stellar collapse 

In the calculation and simulation below, we take DM mass rriy- = lOOGeV, elastic scatter 
cross section 10~^^cm^, inelastic scatter cross section 10~^^cm^ and DM halo energy density 
0.3GeV/cm3 as the fiducial values. 
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3.1 DM distribution and Stellar Evolution 

For our purpose we need to trace the stellar evolution and then the DM distribution inside 
the evolving star. We consider two examples, a 4Mq star as a typical progenitor of white 
dwarfs, and a 15Mq star as a typical progenitor of a neutron star, both at solar metallicity. 
The evolution of both models is followed with the public code MESA [23], using the standard 
numerical methods and input physics as described therein. The 4Mq model has been evolved 
to the formation of a white dwarf, using the mass loss formula for AGB stars by [35], while 
the evolution of the 15Mq model was followed to the formation of the iron core using the 
mass loss formula of [36] for hot stars and of [37] for cool stars. The density profiles of these 
stars at the end of their lives are shown in Fig.l and Fig. 2, respectively. 

The DM distribution, without doubt, changes with the evolution of the stars. The pro- 
cess can be approximated by the adiabatic contraction process. The response of individual 
WIMP to the adiabatic contraction of the star can be traced with the adiabatic invariants 
J and Jj. = 2 /J^™"" \/'^[E — ipi^r)] — J'^/r'^dr, where J is the angular momentum and is 
the radial action. E is the total mechanical energy of the particle and ip{r) is the gravi- 
tational potential at radius r. These invariants are not modified when the star change its 
gravitational potential in the evolution if the particles do not scatter. So we can generate 
the distribution of the WIMPs in every stage of the star evolution with these invariants. The 
detailed description of such a process can be found in [24] and we follow their approach in 
this work. 



On the other hand the DM particles are captured continuously, which lose their energy 
by further scattering and their density profile also changed due to the adiabatic contraction. 
Both effects, the scattering and the adiabatic contraction, are considered in our following 
investigation. 

3.2 Elastic case 

From kinematics we know that a WIMP will lose about ninu/mx of its kinetic energy in each 
elastic scatter. To thermalize a WIMP needs to lose almost all the kinetic energy, i.e.. The 
WIMP will need to scatter roughly m^^/mnu times. Note that the time needed for one scatter 
is — - — . Then the thermalization timescale can be estimated by 

Tther — (o.lj 

PnuCTsV 

where rrinu is the rest mass of the nuclei, p^n is the number density of the nuclei in the 
star, fjs is the elastic scatter cross section and v is the velocity of the WIMPs inside the 
star. We take v to be the average of the escape velocity Ve and the final velocity Vf, where 
Vf == (3A:T/m^)^/^ assuming the dark matter has thermalized with the star and has a 
temperature which is equal to the temperature of the star at its center . in the core of the 
star in our estimation. The central density and the duration of the stars in different burning 
stages are shown in TABLE. 1 in [25], with which it is straightforward to estimate rther and 
then examine whether the thermalization can be established or not. For the stars of 15Mq 
and 4M0 the thermalization can be reached in all these stages. Consequently the DM density 
in the star can be approximated by an isothermal sphere. Before collapse, the thermalized 
WIMPs distribute as [28] 

p(r) = ^ [' ^2{E-ct>{r))eM-E/kT]dE, (3.2) 

2V2 J0(r) 
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Figure 1. The number density profile of the star with an initial mass 15Mq just before its collapse. 



where (j){r) is the gravitational potential before collapse and E is the total mechanical energy 
of the particle. (j){r) and temperature T s evaluated at the end of the star's evolutionary 
track. The baryon density in the core of the star can be taken as a constant, with which 
the integration of Eq.(3.2) give the Eq.(2.1), where the value rth can now be evaluated with 
equation (2.2) and star evolution results. The central T is 3.6 x 10^ K and 1.3 x 10^ K 
for the 15Mq star and 4Mq star, respectively. The density profiles are shown in Fig.l and 
Fig. 2. Then r^h is 10'' cm for the 15Mq and 10^ cm for the AMq star. As the star evolves 
the baryonic profile will adjust. The dark matter profile will also evolve due to adiabatic 
contraction. The distribution function in the energy space is needed, which reads [28] 



p{r,E) 



JdJf{E) 



,_r2[2(E-,^(r))- JVr2]i/2' 



f{E) = (27r)-3/2exp(-i?), 



(3.3) 



(3.4) 



where J is the angular momentum, and Jmax = — (j){r)]/my^ is the maximum angular 

momentum at radius r. 

For white dwarfs, no further collapse happens and the above result is the final profile 
of WIMPs if the annihilation can be ignored. But for neutron stars, the collapse makes the 
situation more complicated. After collapse the gravitational potential has changed a lot, it 
is thus very interesting to see whether a significant fraction of WIMPs is left outside the 
compact object or not. The change of the distribution of WIMPs in the collapse process of 
the medium with a constant density has been discussed in [26]. But in reality the density 



-4- 



le+07 




100000 - 



1000 - 



0.1 



10 - 




le+07 



le+08 



le+09 



le+10 



le+11 



r(cm) 



Figure 2. The number density profile of tlie star witli an initial mass AMq just before the formation 
of a white dwarf. 

profile of tlie dying stars is mucli more complicated (see Fig.l). So we adopt a numerical 
approach. We treat the collapse as a free fall process and calculate the gravity in each time 
step of the collapse, then calculate the corresponding energy change of each dark matter 
particle, until the stellar core has formed a neutron star which has a radius of 12km. After 
the collapse presumably, an explosion happens, which sweeps away material outside of the 
nascent 1.5Mq neutron star. The distribution of WIMPs before collapse is simulated by 
the use of MonteCarlo method with eq.(3.2) and eq.(3.3). We are interested in the WIMPs 
surrounding the compact object. Since the compact object is very dense, any WIMPs passing 
through will scatter with the nuclei so frequently that lose their kinetic energy, and finally 
sink into the center of the compact object. That's why we drop the WIMPs whose semi- 
minor axis is smaller than the radius of the compact object. At the end of the life of the 
4 Mq star, rth ~ 10^ cm is comparable to the typical radius of white dwarfs, so most WIMPs 
are concentrated in the white dwarf and are not of our interest. On the other hand, at the 
end of the life of the 15 Mq star, rth ~ lO'^ cm is only one order larger than the neutron 
star's radius and nearly all the WIMPs falls inside the neutron star in the collapse in our 
simulation. Hence no prominent dark matter mini-halo can be formed. 

3.3 Inelastic case 

The inelastic case [18-20, 29] is different from the elastic case in some aspects. For example, 
an inelastic scatter will not occur if the energy transfer is not large enough. The energy 
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transfer in a scatter is [9] 



AS = -x-nu _ ^-^^^ ^3^^^ 
(m^ + nil, 



''nuj 



where £ is the WIMP's kinetic energy and 9 is the scatter angle. It should be mentioned 
that TTinu ~ 50 GeV since among all nuclei the iron has the greatest contribution to the total 
cross section [18]. In deriving eq.(3.5) the nucleus is assumed to be at rest [10]. Inelastic 
scattering does not happen if £ is below a threshold ifth given by AS = 5, which reads 

m 

fth = (1 + ^)S, (3.6) 
rrinu 

where 6 is the mass gap in the inelastic DM model. For rriy. ~ 100 GeV and 5 ~ 200 keV 
that are consistent with all current experiments [15, 16] we have £th ~ 600keV. 

Generally speaking, the inelastic dark matter should also scatter elastically due to elastic 
coupling [11] or loop corrections if there is no tree level elastic coupling[12]. But the cross 
section of such elastic scattering is highly model-dependent. [11] suggested that the elastic 
cross section can be six orders of magnitude smaller than the inelastic cross section. It is 
shown in [12] that the cross section can be as large as 10~^^cm^ if the force mediator is very 
light while it can also be as small as 10~^'^cm^ if the force mediator is heaver than IGeV. If 
the elastic cross section is large enough the dark matter inside the star is thermalized, for 
which our discussion in last subsection applies and there is no novel phenomena taking place. 
However if the elastic cross section is as small as 10~^^ cm^ the thermalization can not be 
reached in our star evolution model. Such a scenario is the focus of the following discussion. 
With our Eq.(3.1) and Table 1 of [25] it is straightforward to show that the thermalization 
of dark matter particles is established most easily in the helium burning stage. Therefore we 
have a threshold of the elastic cross section below which the dark matter particles cannot be 
thermalized in all burning phases. For our 15M0 model the threshold is ~ 10~^^ cm^ while 
for our 4Mq model the threshold is ~ 10~^^ cm^. As a result, if the elastic cross section is 
very small, say, about 10~^° cm^ as mentioned above, thermalization cannot be reached. We 
assume it is the case in our following discussion. 

Due to the inelasticity, WIMPs do not thermalize with the star and do not have a 
thermal distribution as in elastic case. A WIMP will lose most of its kinetic energy after 
sufficient scattering and then it has a residual kinetic energy below iSth, which depends mainly 
on the kinematics of the last scattering. So the DM density profile can not be treated as an 
isothermal sphere any longer. A simple assumption is that the kinetic energy of the WIMPs 
after enough scattering has a flat distribution ranging from zero to 6. If the possibility of 
the scattering at different radius is known, the distribution of total energy E can be derived 
by adding the kinetic energy and gravitational potential. But note that the possibility of 
scattering at different radius depends on the particle's orbits, which will change significantly 
in former scatters. So in this work we carry out a Monte-Carlo simulation to trace the 
particles' fate inside the star. We calculate their orbits after each scattering and then simulate 
the next scattering with the kinematics described above. The evolution of stars of 15Mq and 
4:Mq has been divided into tens of intervals. We calculate the DM capture rate and the 
resulting DM density profile in each interval and also take into account possible scatter and 
adiabatic contraction in the later phases. The rest of the simulation, i.e., the DM density 
change due to the collapse, is similar to the elastic case. Please see figure 3 and figure 4 for 
the results for a neutron star and a white dwarf, respectively. In the case of a neutron star 
there are a large amount of WIMPs outside the compact object, i.e., a DM mini-halo with 
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Figure 3. The density distribution of the DM mini- halo surrounding the neutron star: the inelastic 
DM model. 

a radius of ~ 2 x 10^ cm is formed. In the case of white dwarf, the radius of the compact 
object is only a few times smaller than the size of the DM mini-halo. Hence no prominent 
mini-halo forms surrounding the white dwarfs. 

Please note that in the simulation we use the fiducial DM halo density 0.3GeV/cm'^. In 
some regions of the galaxy, such as the galaxy center and some DM clumps, the DM density 
can be much higher. If the lifetime of the star is much smaller than the equilibrium time 
Teq, the annihilation of DM inside the star is not important compared with the capture. So 
the DM density of the mini-halo is proportional to the capture rate and then the local DM 
density. If before collapse the equilibrium has already been established, the DM density in the 
mini-halo should be proportional to the equilibrium density. In such a case the capture rate C 
is equal to the annihilation rate CaN'^, which means that the DM density is proportional to 
the square root of the local DM density. Therefore in both cases the denser the environment 
DM, the higher the DM density of the mini-halo. 

3.4 General relativity correction 

Near the compact objects, in particular the neutron star, the gravitation is so strong that 
general relativity may be not ignorable any longer. In the simulation we find that the captured 
particles all have a very large eccentricity, so most orbits lie far from the central star, where 
the general relativity effect can be ignored. The general relativity effect we are interested 
in is that the particles moving around the central star would have no stable orbits if their 
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Figure 4. The density distribution of the DM mini-halo surrounding the white dwarf: the inelastic 
DM model 



angular momentum is small enough. As shown in [22], if the condition 

J<2V3^^ 



(3.7) 



is met the particles will be captured by the central object in a Schwarzschild metric, where 
M is the mass of the compact object and c is the speed of light. Other less important general 
relativity effects such as the distortion of orbit shape and time delay have been ignored. 

3.5 Recoil velocity 

When the star collapses a huge amount of energy is released. The energy release is not 
always isotropic and the formed compact object may get a considerable recoil velocity. For 
the neutron star the typical velocity is about 300 km s^^ and for white dwarf it is much 
smaller. Since the collapse is transient, the recoil can be taken as a 'kick' of the particles 
in the compact object's frame, i.e., all particles captured by the star get a velocity at the 
same direction and their kinetic energy change. When the collapse happens the gravity is 
so strong that the escape velocity is of the order lO^km/s , so a small change of the kinetic 
energy can be safely neglected in our simulation. 

4 Time evolution of the density of the DM mini-halo 

As mentioned above, only in the case of inelastic scattering for a neutron star does an 
extended mini-halo form. Below we focus on such a specific scenario. A WIMP in the DM 
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mini-halo will maintain its orbit except annihilating with the other WIMP or scattering with 
baryons. Since we consider only the WIMPs whose orbits are totally outside the neutron star, 
the scattering with baryon can be neglected. The time evolution of the density of WIMPs is 
governed by 

rl n 

- < aAV > p . (4.1) 



dt 

The solution is simply 

P=—, , (4.2) 

1/P0+ <aAV>t' ^ ^ 

where po is the number density at time 0, i.e., the moment right after the collapse. The 
supernova remnant surrounding neutron star persists for about ~ 10^ years, after which 
the remnant material becomes transparent and the products of WIMPs annihilation become 
detectable. For t ~ 10'^ years and < a^v >~ 10 cm'^s ^, number density p won't change 
significantly unless it is larger than 10^^ cm~^. For the cosmological timescale t ^ 10^^ 
years, the DM number density evolve due to annihilation can be as high as ~ 10^ cm 
From fig. 3 we find that the DM number density in the nascent mini halo is much larger 
than ~ 10^ cm ^. As a result of the annihilation, the number density drops with time as 
p~ 10^°(t/10^ year)-^ cm'^ for t > 10^ years. 



5 Flux of annihilation products 

Now we can estimate the flux of the annihilation products of WIMPs in the mini-halo sur- 
rounding the neutron stars. The annihilation rate is 

iVanni = j < ^av > p^{r)r'^dr. (5.1) 

In the inelastic DM model, as shown in flg.3, the WIMPs are centered in a small region 
r ~ 2 X lO^cm with the number density p{r) ~ 10^^ cm~^ shortly after the collapse. The 
current number density is in order of 10^ cm~^ and the annihilation rate is ~ 10^^ s~^ if the 
compact objects are old enough. Let's estimate the contribution of the Galactic DM halo. 
The local DM energy density is 0.3 GeV cm~^, for ~ 100 GeV the annihilation rate of 
Galactic DM halo particles in ~ 1 kpc^ space is about 10^^ per second. The number densities 
of local neutron stars is ~ 2 x 10^ kpc"'^ [22]. So for inelastic WIMPs the annihilation rate 
of these mini-halos in 1 kpc'^ is about 10^^ s~^, i.e.. The amount of DM annihilating in the 
mini-halo is small compared to the Milky Way Halo. 

It is interesting to investigate whether the annihilation of inelastic DM particles in some 
nascent mini-halos can produce enough products giving rise to detectable signal. To account 
for the cosmic electron/positron data at an energy ~ 100 GeV, the annihilation rate should 
be ~ 10^^ — 10'^^ s~^ kpc~^. Therefore no observable signal is expected unless one nascent 
mini-halo is within a distance ~ a few 0.1 pc to the Earth. 



6 Conclusion 



In this work we have investigated the formation, evolution and the annihilation of the particles 
in dark matter mini-halo surrounding the compact object. In our model, the dark matter 
particles are captured by the stellar matter and then distributed in the core of the star. At 
the end of its life the star uses up the fuel and a compact object, either a white dwarf, a 
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neutron star or a black hole, is formed. The dark matter particles will redistribute around 
the newly formed compact object. In the elastic dark matter model and if the elastic cross 
section is large enough in the inelastic dark matter model, the dark matter particles will be 
thermalized and no prominent dark matter mini-halo can be formed. In the inelastic dark 
matter model, however, if the elastic cross section is very small (for example < 10"^'' cm^) 
an extended dense dark matter mini-halo surrounding neutron star may be formed. Though 
such mini-halos may be very common in the Galaxy, the electron/positron flux resulting in 
the annihilation of these dark matter particles is negligible compared with the observational 
data unless a nascent mini-halo is within ~ a few 0.1 pc to the Earth. 
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